o 

(N 



Neutrino Dark Energy With More Than One Neutrino Species 

Ole Eggers BjaeldeQ 
Department of Physics and Astronomy, University of Aarhus, 
Ny Munkegade, DK-8000 Aarhus C, Denmark and 
Center for Cosmology and Particle Physics, Department of Physics, New York University, New York, NY 10003 

Steen Hannestad 

Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, DK-8000 Aarhus C, Denmark 

(Dated: May 17, 2010) 

The mass varying neutrino scenario is a model that successfully explains the origin of dark energy 
while at the same time solves the coincidence problem. The model is, however, heavily constrained 
by its stability towards the formation of neutrino bound states when the neutrinos become nonrela- 
tivistic. We discuss these constraints and find that natural, adiabatic, stable models with the right 
amount of dark energy today do not exist. Secondly, we explain why using the lightest neutrino, 
which is still relativistic, as an explanation for dark energy does not work because of a feedback 
mechanism from the heavier neutrinos. 
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I. INTRODUCTION 



II. MASS VARYING NEUTRINOS 
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Precision observations of the cosmic microwave back- 
ground [H-[3] , the large scale structure of galaxies [3| , and 
distant type la supernovae have led to a new stan- 
dard model of cosmology in which the energy density is 
dominated by dark energy with negative pressure, lead- 
ing to an accelerated expansion of the universe. 

A very interesting proposal is the so-called mass vary- 
ing neutrino (MaVaN) model in which a light 
scalar field couples to neutrinos, see also (i~2rl5l| for back- 
ground. Because of the coupling, the mass of the scalar 
field does not have to be as small as the Hubble scale but 
can be larger, while the model still accomplishes late-time 
acceleration. 

In this paper we discuss the different criteria that need 
to be fulfilled in order to have a stable, adiabatic mass 
varying neutrino model with the correct cosmology to- 
day - i.e. Ode = 0.7, Qm = 0.3 and w ~ —1. Our aim 
is to show that it is only possible as long as the scalar 
field potential resembles a cosmological constant - and 
hence the model loses its prediction of the current neu- 
trino mass and fails to solve the coincidence problem. 

In addition, we discuss the suggestion that the lightest 
neutrino, which can be relativistic today, may be respon- 
sible for dark energy. We find that there is evidence that 
the relativistic neutrino will feel an instability towards 
the formation of neutrino nuggets. 

In the next section we briefly review the formalism 
needed to study mass varying neutrinos and in Sec. IIIII 
we discuss the different criteria in the MaVaN framework. 
In Sec. lIVI we discuss MaVaNs with a relativistic neutrino 
and in Sec. [V] we conclude. 



The idea in the mass varying neutrino scenario [9H11I] 
is to introduce a coupling between neutrinos and a light 
scalar field and to identify the coupled fluid with dark en- 
ergy. In this scenario the neutrino mass m„ is generated 
from the vacuum expectation value (VEV) of the scalar 
field. Thus at scale factor a the pressure P v {m v {(j)),a) 
and energy density p v {rn v {^) , a) of the uniform neutrino 
background contribute to the effective potential V(<p, a) 
of the scalar field in the following way 



V{cj>) = V^) + (p v - W v ) 



(1) 



where V^,(<p) denotes the fundamental scalar potential. 
The energy density and pressure of the scalar field are 
given by the usual expressions, 



(2) 



Defining w = Pde/pde to be the equation of state 
of the coupled dark energy fluid, where Pde = P v + P<b 
denotes its pressure and pde = Pv + P<t> its energy density, 
the requirement of energy conservation gives, 



Pde + 3Hp DE (l + w) = 0. 



(3) 



Here H = - and we use dots to refer to the derivative 

a 

with respect to conformal time. Taking Eq. ([3]) into ac- 
count, one arrives at a modified Klein-Gordon equation 
describing the evolution of <fi, 
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Here and in the following primes denote derivatives 
with respect to (' = d/dcj)) and (3 = feg2s is the 
coupling between the scalar field and the neutrinos. 

In the nonrelativistic case P v ~ 0, such that Eq. (p} 
takes the form 



^CDM^CDM 

nient form 



fthSh. This can be recast in a more conve- 



-f2„ < J7m, where we have ne- 



glected the effect of baryons compared to cold dark mat- 
ter and we have assumed the density contrasts of roughly 
the same order. 



V = p„ + Vj, 



(5) 



Assuming the curvature scale of the potential and thus 
the mass of the scalar field m$ to be much larger than 
the expansion rate of the Universe, 



III. MODEL REQUIREMENTS 

In summary we have the following different criteria for 
an adiabatic MaVaN model. 



V" 



(6) 



the adiabatic solution to the equation of motion of the 
scalar field in Eq. d3| applies As a consequence, 

the scalar field instantaneously tracks the minimum of 
its effective potential V 



V = p' v + Vi=(3 Pu + Vi=0 



(7) 



MaVaN models has a possibility of becoming unstable 
on sub-Hubble scales m^ 1 < a/k < H^ 1 in the nonrela- 
tivistic regime of the neutrinos, where the perturbations 
5p v evolve adiabatically. 

In Ref. it is shown that the equation of motion 
for the neutrino density contrast [62j in the regime 
mT 1 < a/k < H^ 1 can be written as 



= -H* 



CDM»CDM 



rt h s, 



(8) 



where 



G e ff 

and 



G 1 



2(3 2 M 2 



pi 



a 2 Pi 
3ff 2 M 2 ' 



(9) 
(10) 



Since neutrinos not only interact through gravity, but 
also through the force mediated by the scalar field, they 
feel an effective Newton's constant G e g as defined in 
Eq. © . The force depends upon the MaVaN model spe- 
cific functions (3 and V$ and takes values between G and 
G(l + 2/3 2 A/ 2 [) on very large and small length scales, re- 
spectively. 

In certain cases of strong coupling neutrinos suffer an 
instability towards clumping in which case they stop be- 
having as dark energy [2l[ . In Ref. [13] a criterion for the 

2^ Ml, 



stability was developed. 
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0, U 5 V < 



1) The model needs to satisfy current observations 
i.e. w *~ — 1. This is easily fulfilled by demand- 
ing <fi r~j since we know the neutrino contribution 
to both pressure and energy density is smaller than 
the scalar field contribution. 

2) We want a neutrino mass at present that satisfies 
the current neutrino mass bounds from cosmology, 
that is m„ < leV. Using this, we get a maximal 
value for the neutrino density parameter f2„ < 0.02 

aim. 

3) The model must produce the right cosmology, i.e. 
VL DE ~ 0.70. This gives us roughly ^ - ^ > 35. 

4) In order to have a prediction for the current neu- 
trino mass, we are looking for an adiabatic model 
that continuously tracks the minimum of the effec- 
tive potential, i.e. V = VI + fip v = 0. In addition, 
to avoid severe fine-tuning, we demand m^ 1 <C 
H- 1 [HEH. Hence V$ + (/3 2 +P')p v > H 2 , using 
a constant coupling gives us VI' + 1 p v 3> H 2 . 

5) We want our model to be stable towards the for- 
mation of neutrino nuggets at present. This gives 



the requirement 



Q. u < 51m- Since we 



are in the adiabatic regime, we are looking on scales 
where k/a ~ yjvf + Pvfi' ■ Combined with the fact 
that Qm ~ 0.30, this criterion simply reduces to 



P 2 M 2 Y < 15. 
Combining the criteria above and using H 2 



Ptutul 

"M 5 " 



-jjpr , we can see that in the case of a constant coupling we 
pi 



1 or roughly 



72— r 



which reduces to 



15 

35 



V" M 2 



> 1. 
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We propose potentials that resemble those originally 
proposed in Ref. (Tlj . and in order to be able to explain 
the coincidence problem we want a potential as simple as 
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possible without any type of cosmological constant [63]. 
However all single field potentials previously studied in 
this context do not fulfill the constraint Eq. (fTTj) . 

Hence we need a sort of hybrid model as suggested 
by (26|. In this type of model we have effectively two 
minima. One is the true minimum which has zero energy 
density and the other is the false minimum in which our 
current universe sits. The offset between the two minima 
is then what is referred to as the dark energy density. The 
dynamics of the scalar field is given by the movement of 
the false minimum with the dilution of the universe. One 
very nice feature of this model is that it is implemented 
in supersymmetry and hence explains the stability of a 
very small scalar field mass. However, since we have an 
offset it can be interpreted as an effective cosmological 
constant and as such the hybrid model fails to address 
the coincidence problem. In addition, to keep the model 
stable it is suggested that the lightest neutrino is still 
rclativistic and that we associate dark energy with this. 
We address this issue below. 



IV. MAVAN MODEL WITH A RELATIVISTIC 
NEUTRINO 

We aim at keeping the model as simple as possible, so 
we propose the existence of just one scalar field with a 
democratic coupling to all light neutrino species (64j . In 
that case it is assumed that the lightest neutrino does not 
clump because of pressure, while the two heavier neutri- 
nos become unstable. In the following we will argue that 
this is not possible since a feedback from the growth of 
the heavy neutrino perturbations will cause the relativis- 
tic neutrino perturbation to grow as well. The problem 
in this scenario can be neatly illustrated by the following 
equation[22j 



Spu = ^ / q 2 dqdnef (q)^ + 5^{p v - 3P„). 



(12) 



The equation explains the growth of the neutrino per- 
turbation and applies to both the interaction between a 
relativistic neutrino and a scalar field as well as that of 
a nonrelativistic neutrino and a scalar field. 

Let us consider a system consisting of three neutrinos, 
two heavy nonrelativistic neutrinos (H) and one light rel- 
ativistic neutrino (L) both interacting with a scalar field. 
We only have one scalar field so the nonrelativistic neu- 
trinos will inevitably become unstable towards clumping 
and 8p v {H) grows to very large values. For an interac- 
tion with a nonrelativistic neutrino, the average scalar 
field perturbation can be written as [27] 



(Vl' + P^') 



fc2 



(13) 



As we can see the scalar field perturbation is effectively 
proportional to the neutrino perturbation. Note that in 



reality there will be extra terms in Eq. ([13]) since we also 
have a relativistic neutrino, however, these terms will be 
sub-dominant compared to the terms from the heavier 
neutrinos. So 50 will grow as a result. 

Turning our attention to Sp u (L), we consider the sec- 
ond term in Eq. (j L2|) . 8<f) is growing uninhibited to large 
values, the coupling ft is for simplicity assumed of or- 
der unity (Planck units) [65]. (p v — 3-P„) is a suppression 
factor of the order m/E - this factor will act to delay 
the growth of 8p u (L). However, since 8<j> will continue 
its growth, the inevitable conclusion is that 5p u (L) will 
eventually start to grow. Hence, there exists a type of 
feedback mechanism between the neutrinos. Note that 
in this treatment we only considered the interaction be- 
tween neutrinos of the same type via exchange of the 
scalar field. If we include couplings between different 
neutrinos, the feedback mechanism will be even stronger. 
One could of course argue that we are exactly living in 
a transition regime when 8p„(L) has still not turned un- 
stable. However, that would require serious fine-tuning. 

We also investigate a MaVaN model with a relativistic 
neutrino numerically using the publicly available CMB- 
FAST code (54|, where we implement the scenario with 
two heavy neutrinos with assumed masses m v {H) = 
0.312 eV and one species of relativistic neutrinos with 
assumed mass m v (L) = 0.0001 eV. For the scalar field 
potential we use a Coleman- Weinberg [55j type scalar field 
potential similar to the one presented in Refs. [TJ H3 



= V log (1 + nq 



(14) 



where re and Vq are dimensionful constants, which are 
fixed by demanding the appropriate amount of dark en- 
ergy in the universe today as well as solving for the min- 
imum of the effective potential. For the mass term of the 
three neutrinos we choose the same mass term as in 127] 



mo 



(15) 



which results in a coupling /? = — i [66[ . 

The code calculates the background energy density and 
pressure of different species from standard integration 
methods while at the same time it solves for the minimum 
of the effective potential of the neutrino-scalar field fluid 
which is 



V{4>) = W) + (PAH) - 3Pu(H)) + ( Pv (L) - 3P„(L)), 

(16) 

where we remember that there are now two heavy (H) 
neutrinos and one light (L) neutrino. The expression 
above can be recast in a more convenient form as 12611 



V{(j>) = V^{cj>)+m v {H)n v {H) 
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(17) 
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where we introduced the neutrino number density n v 
and neutrino temperature T v . At the same time as track- 
ing the background we calculate the perturbations to first 
order such as 



pv = Pv + 5p v . 



(18) 



This is done for both heavy and light neutrinos using 
Eq. (fT2|) . The results are shown in Fig. [T] where we 
plot the evolution of the neutrino perturbation and CDM 
perturbation from a redshift of z + 1 = 11 until today. 

What happens in Fig. Q] is that the perturbations of 
the heavy neutrinos, which have already turned nonrela- 
tivistic before the range of this plot, are already growing 
albeit at a slow pace because of the growth-slowing ef- 
fect of CDM - see Eq. ©. The perturbation of the light 
neutrino is oscillating (only visible if zoomed in) in this 
epoch, but also growing slightly as a result of the growth 
of the coupling f3. The growth-slowing effect of CDM is, 
however, only temporary and around a reshift of z+1 ~ 2 
the heavy neutrino perturbation can no longer resist the 
temptation towards unstable growth. The heavy neu- 
trino perturbation explodes and because of the increased 
value of /3, the light neutrino perturbation starts explod- 
ing as well. We do see that eventually the CDM per- 
turbation also starts exploding and at this moment the 
linear code is no longer valid. However at this moment 
both neutrino perturbations are already in the nonlin- 
ear regime and will have started forming neutrino bound 
states. In this regime the MaVaN fluid will no longer 
act as dark energy and the model breaks down. In con- 
clusion, as a result of a feedback mechanism, the blowing 
behaviour of the nonrelativistic neutrino density contrast 
causes the relativistic neutrinos to start clumping as well. 
Hence the neutrino scalar-field fluid will start acting as a 
cold dark matter component (clustering neutrinos) and 
hence cannot be attributed to dark energy. 

The reason that the relativistic neutrino is able to 
clump is that it will acquire an effective mass, thus it can- 
not be regarded as a relativistic particle. Unfortunately, 
we cannot use conventional bounds to constrain this ef- 
fect since the evolution of the neutrino perturbations be- 
come nonlinear, i.e. the whole system of equations we are 
solving, starting with the modified Klein-Gordon equa- 
tion breaks down. This has the effect that all current 
bounds are no longer valid, as these are established in 
the linear regime. 
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FIG. 1: Density contrasts plotted as a function of redshift for 
a system consisting of one light and two heavy MaVaN neu- 
trinos each interacting with the same scalar field. The scale 
is k = 0.1 Mpc -1 and we choose the current neutrino masses 
m v (L) = 0.0001 eV and m v (H) = 0.312 eV (Note that the 
choice of current neutrino masses does not affect the result 
qualitatively). We have chosen k = 1 x 10 18 in CMBFAST 
units of (Af PL /Mpc)- 1/2 and V = 1.11 x 10~ 9 (M PL /Mpc) 2 
to fix the cosmology. The solid black line is CDM-density 
contrast defined as <5pcDM , the dotted blue line is the heavy 

PCDM ' J 

neutrino density contrast iEidSl^ anc j the dashed red line is 



the light neutrino density contrast 



The heavy neutri- 



nos grow moderately until the coupling becomes large enough 
for the instabilities to set in. The light neutrino is still rela- 
tivistic, and its density contrast oscillates as acoustic waves 
while growing slightly. However, due to a feedback mech- 
anism, the relativistic neutrino density contrast tracks that 
of the nonrelativistic neutrino especially as the growth of the 
heavy neutrino perturbation becomes unstable; i.e. both neu- 
trino species will clump. Note that the CDM perturbations 
also blow up at late times. This is an effect of the system of 
differential equations breaking down as all parameters go to 
infinity. 



potential resembles that of a cosmological constant. In 
this case we simply move the fine-tuning problem to that 
of explaining the small fraction in the power-law - which 
unfortunately is also rather unnatural. 



V. CONCLUSION 

Single scalar field models can be used to explain late- 
time acceleration in the MaVaN scenario. However, in 
general using these potentials leads to instabilities to- 
wards neutrino bound states unless certain criteria are 
relaxed. In order to obtain a stable or metastable model 
we direct the reader to Ref. [13] in which one such model 
is presented. The crucial thing is that the scalar field 



Accordingly it has been suggested to include an ex- 
tra scalar field in the treatment. This has some very nice 
features and is easily capable of obtaining late-time accel- 
eration as well as fluE = 0.7 today. However, one draw- 
back is the need for the lightest neutrino to be relativistic 
today. As was explained above the feedback mechanism 
will eventually cause the perturbations of the relativistic 
neutrino to start growing. As a result, the MaVaN fluid 
will cease to act as dark energy. 
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